Scalable register initialization for quantum computing in an optical lattice 
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The Mott insulator state created by loading an atomic Bose-Einstein condensate (BEC) into an optical lattice 
may be used as a means to prepare a register of atomic qubits in a quantum computer. Such architecture requires 
a lattice commensurately filled with atoms, which corresponds to the insulator state only in the limit of zero 
inter-well tunneling. We show that a lattice with spatial inhomogeneity created by a quadratic magnetic trapping 
potential can be used to isolate a subspace in the center which is impervious to hole-hoping. Components of 
the wavefunction with more than one atom in any well can be projected out by selective measurement on a 
molecular photo-associative transition. Maintaining the molecular coupling induces a quantum Zeno effect that 
can sustain a commensurately filled register for the duration of a quantum computation. 
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In the past decade, tremendous progress has been made 
in the field of atomic physics toward the creation of macro- 
scopic quantum states such as Bose-Einstein condensates on 
one hand id and precise control of single and few coupled 
atoms on the other, for example in optical lattices [2]. Several 
years ago it was proposed 1 3] to marry these advances by load- 
ing an optical lattice from an atomic BEC. If one begins with 
a superfluid-like BEC and adiabatically turns on an optical 
lattice potential, the system will undego a phase transition to 
the Mott insulator (MI) state (characterized by the same num- 
ber of atoms in each lattice well) when the intra-well interac- 
tion energy is much greater than the inter-well tunneling. 1 4). 
Recent experiments 1 5 ] have demonstrated the Mott insulator 
phase transition in a magnetically confined three dimensional 
optical lattice, with an average filling factor of two atoms per 
well. In this paper we will show how the use of an quadratic 
trapping potential and selective measurement of atomic pairs 
allows for the MI transition to become a robust mechanism for 
quantum register initialization. 

It was recognized early on that the MI transition might be 
an efficient way to initialize a register of atomic qubits in an 
optical lattice for use in quantum information processing. A 
key advantage of loading from a BEC is the availability of an 
initially high phase space density which can be frozen to the 
MI state with atoms occupying every lattice site. When the 
lattice is loaded such that only the lowest vibrational state of 
each lattice well is occupied, the system is well described by 
the Bose-Hubbard Hamiltonian: 

Hbh =Y, e U) n j- J ( a j a j+i + a )+\ a i) + W n i( n j - 1) (!) 
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Here aj are the bosonic annihilation operators and nj = a -aj 
the number operators for an atom in the lowest vibrational 
state of lattice well j. The energy offset at each lattice site 
is which models a continuously varying external po- 
tential. The energies J and U are the tunneling and on- 
site interaction energies respectively. In the tight-binding 
model, the nearest neighbor tunneling energy J is defined 
as one fourth the band width of the lowest occupied band. 
For tunneling through a potential barrier given by V(x) = 
Vcos 2 (fcic), the tunneling rate is closely approximated by [6] 



J/h = 4/(y/nh)E R (V /£ s ) 3 / 4 e~ 2 V^7£«, where the recoil en- 
ergy is Eg = (hk) 2 /2m (m=atomic mass). The on site inter- 
action is a result of the ground state collisions described by 
the s-wave scattering length a s between two atoms each in the 

motional state 0(x) and is given by U = ^jjp- J dx\$(x)\ 4 . 

For the homogeneous system (e(j) = Vj) of fixed extent, 
the behavior of the system is uniquely described by the ra- 
tio U /J which for trapping in an optical lattice decreases ex- 
ponentially with the trap depth V. While our results are ap- 
plicable to higher dimensions, henceforth we assume a three 
dimensional lattice with tight transverse confinement and tun- 
neling dynamics along one dimension only. For the homo- 
geneous system, only commensurate fillings give rise to a MI 
transition. However, one should note that an adiabatic transfer 
mechanism between two sublevels of each atom may be used 
to fix nonuniform filling |7]. A caveat is that even with unit 
filling, the MI state still carries small but non-zero number 
fluctuations which provide a small residual coherence across 
the system that scales as the number of trapped atoms 
Applying first order perturbation theory in Hbh, the ground 
state for N atoms in M wells in one dimension with N = M is 
approximately 

\V g )=a(\T)+2VNJ/U\S)), (2) 

where the normalization constant is oc= (1 + AN(J /U) 2 )~ x l 2 . 
Here the unit filled target state is \T) = IlJ=i fl J|0) and the 
symmetrized state, assuming periodic boundary conditions 
(j + M = j), is \S) = l/V2NZjl 1 (a f J+1 aj + a)a j+1 )\T). The 
energy of the ground state is approximately E g = —ANJ 2 /U. 
In general, a homogeneous lattice will have a probability of 
approximately l/M of being commensurately filled and thus 
does not provide a robust mechanism for initializing a quan- 
tum computer. 

We therefore propose to use an inhomogeneous lattice with 
open boundaries created by a weak quadratic magnetic trap 
that acts to collect atoms near the center of the trap and leaves 
empty wells (holes) at the edges. For our analysis we assume 
a one dimensional optical lattice, with N < M | 10], in the pres- 
ence of a weak magnetic trap with oscillation frequency (Oj. 
The characteristic trap energy scale 8 = m/2{%/k) 2 (£>\ is de- 
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fined so that e(j) — 5j 2 . We stipulate that the on site interac- 
tion energy be larger than the trapping energy of the most ex- 
ternally trapped atom, 01U > e( (N — 1) /2), in order to inhibit 
multiple atom occupation in any well. The register is defined 
by a physical subspace comprising a number of wells n <N 
in the center region of the trap. The barrier space flanking 
will act to suppress percolation of holes from the edges to the 
center. The estimated probability for holes in ^ due to tun- 
neling through the barrier is p h » U%^lf) /2 (^/ 8 ( 2 J + 1 ) ) 2 = 
(7/25) A '-" +2 (r[n/2]/r[^/2 + l]) 2 , which is negligible pro- 
vided the barrier region is sufficiently large and J/nb < 1. 

When expanded in the Fock state basis, the ground state 
of the register has amplitude in those states with holes neigh- 
boring atomic pairs, analogous to the homogeneous case. We 
describe a protocol which projects out these components by a 
null result from selective measurement of atomic pairs within 
any lattice site. This measurement detects population on an 
excited molecular state and can be made with high efficiency. 
Once the unit filled state is reached with high confidence, con- 
tinuing the measurement will maintain this state by virtue of 
the quantum Zeno effect. 

The measurement will map the register from the ground 
state I'Fg) to the unit filled target state which is not an eigen- 
state of Hbh- To describe the dynamics in the register during 
the measurement, it is convenient to use the following incom- 
plete basis over : 



(n-i)/2 
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For each j the states \Sj) are distinguished by the two ener- 
getically distinct orientations of an atomic pair and its neigh- 
boring hole with energies, E(Sj) = U(l =F jj(2./' — 1)). The 
target state \T) defines the zero of energy. The \T) and the 
{Sj 1 ) states are coupled to first order in Hbh and they span the 
reduced state in Hi of the ground state of the total system. 

In the limit of large n, the dynamics of the register is re- 
stricted to the basis of Eq. [3] This argument is understood 
by comparison to the dynamics in the homogeneous sys- 
tem. In the latter, the state with the largest coupling from 
the target state is the symmetrized \S) with coupling ma- 
trix element (T\Hbh\S) = —2^/nJ. The state \S) itself cou- 
ples to a symmetrized state \S') of all Fock states with a one 
site separation between the atomic pair and the hole: \S') = 
1 / v / 2nLj( fl ) fl ;+2 + a )+2 a j)\T)- The coupling between these 
states is (S'\Hbh\S) = —37. The dynamics on time scales 
t < l/J are therefore constrained to the subspace 
when n> 1. In the inhomogeneous case, the degeneracy 
is absent between states with neighboring pairs and holes, 
{15^)} and states where pairs and holes are separated. In 
this off resonant situation, the coupling to states outside the 
restricted subspace can only be smaller than in the homoge- 
neous case. 

For the measurement, we choose a catalysis laser that is 
on resonance from the ground state of two atoms in a sin- 
gle well to a bound state v of a dipole-dipole coupled molec- 
ular S + P state. The bound state is chosen such that the 
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FIG. 1 : Schematic of the relevant couplings in the problem. The unit 
filled state \T) describing a target quantum register and the states 
1 5^) having one doubly occupied lattice site and a neighboring hole 
are coupled to first order in Hbh ■ A catalysis laser couples states | ) 

to the excited states ) having a bound molecule at the doubly oc- 
cupied site. The bound states quickly decay and give the possibility 
of monitoring population in the "faulty" register states \S^). 



catalysis laser is far off resonance from other bound states 
and repulsive potentials, see Fig[2 For our many body 
system, we adopt the set of many body states {\M^) — 

(l/V2)fct +(1=Fl)/2 a2 +(lTl)/2 |S^)}, where is the creation 
operator for a molecule in the bound state v. These states de- 
scribe n — 2 atoms trapped in the lattice and a single molecule 
at site j + (1 =p l)/2, with dipole-dipole coupling energy 
{M^\H ( [ c i\M j ) = h(Oy. The free atomic Hamiltonian for n 
atoms is Ha = Y.j^egle^iej] where \ej)(\gj)) denotes the 
excited(ground) state for an atom at site j. Hereafter, we work 
in units with the numerical value of h equal to 1 meaning en- 
ergies are understood as being in units of inverse time. The 
"bare" energy Hamiltonian Ho including coupling in the re- 
stricted basis of Hbh is then: 



H = H A +H dd + H BH 

= CO^L,- \ ej ){ej\ +^ ± E(Sf)\Sf)(Sf 
+ ((Ov+E(Sf)-U)\Mf){Mf\ 
-V2JLjA\Sf)(T\ + \T}(Sf\), 



(4) 



Under the atom laser interaction, Hal, the ground and ex- 
cited state of each atom is coupled as is each many body state 
|M*) and its corresponding ground state 15^). In the rotating 
wave approximation, the interaction is: 



Hal - %-Lj(e-*»\ej)<gj\+h.c.) 

+^ ] Ae-'^'\Mf){Sj\+h.c.) 



(5) 



where Qa{m) are the atomic (molecular) Rabi frequencies, re- 
lated by Q.m = VTv^a, where F v is the Franck-Condon factor 
equaling the spatial overlap between the relative coordinate 
wavefunction describing two ground electronic state atoms 
trapped in a single lattice well and the molecular bound state 
V. For bound states of interest, such as the long range bound 
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states of the 0~ (P3/2) potential lUTIl . the catalysis detun- 
ing from atomic resonance, A = co^ — (O eg , is several thou- 
sands linewidths meaning the atomic saturation is low sa — 
(£2^/2 A 2 ) <C 1. In this case, the excited atomic states can be 
eliminated and the total Hamiltonian Hq +Hal in the rotating 
frame is, 



Hi = Lj.±m\+E(Sf))\Sf)(Sf\) 
+(\Vc\+E(Sf)-U)\Mf)(Mf\ 
-V2J(\Sf)(T\ + \T)(sf\) 
+ ^(\Mf)(Sf\ + \Sf)(Mf\)). 



(6) 



Here, the differential light shift seen by the \T) state and the 
\Sj) and \M^) states is equal to twice the single atom light 
shift: V c = Asa- Any population in the bound molecular states 
will decay at a rate Jm ~ 2T, where F is the single atom decay 
rate. Because the excited state linewidth is greater than the 
lattice depth, the decay products will typically not be trapped 
and the system evolves according to a trace non-preserving 
master equation: 

p = -i[H h p] -y M /2^(|M±)(M±|p + p|M ; ± )(M±|). (7) 

j,± 

We have ignored spontaneous emission due to decay from the 
single atom excited states at a rate saT <§; F per atom. For time 
scales 1 /ym <f < 1 /(U + \V C \), the excited state coherences 
can be solved for. To second order in U/Jm, \V c ■{/ju an d first 
order in J/Jm they are 
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(8) 



Inserting these expressions back into the equations for the dy- 
namics in the ground state we have 
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These equations describe the Bose-Hubbard coupled states 
with a decay in population of each state with an atomic pair 
at a rate 2k = Q^Jm / (4(t/ 2 + (Ym/2) 2 ), and a decay of the 
coherences between each of these states and the \T) state at 
approximately the rate K. 

This type of evolution characterized by measurement in- 
duced phase damping was studied extensively by Gagen and 
Milburn IU2I1 . We now show that our system can satisfy the 
conditions for this effect and in particular can be driven to the 
I T) state by monitoring the environment for a signature of de- 
cay from the molecular bound state. 

For the inhomogeneous system, the state \T) couples to 2n 
distinguishable states \Sj). However, we can define an effec- 
tive Rabi frequency between the subspace \T) and the sub- 
space spanned by {|5^)}. This frequency is close to the cou- 
pling matrix element between the state \T) and the state \S) 



in the homogeneous system, namely 2^/nJ. The coherences 
between the two subspaces decay at a rate K, and the popula- 
tion in the subspace {l^j 1 )} decays at a rate 2k. The "good" 
measurement regime as derived in 1 12] is then: 



-« 1 <y^7- 

1m WnJ 



(10) 



The left side inequality ensures that the excited states \Mj) 
are weakly populated (equivalent to the condition for adiabatic 
elimination of these states). The right side inequality ensures 
that measurement is sufficiently strong to damp coherences on 
the time scale that they develop due to tunneling. 

The limiting quantity that determines the decay rate of the 
weakly saturated molecular states and hence the measurement 
strength is the Franck-Condon factor F v . It is calculated for 
bound-bound transitions using the reflection approximation 
of Julienne [14] where we solve for the ground state rela- 
tive coordinate wavefunction for two atoms in a lattice well 
using a pseudo potential appropriate for &1 Rb. We choose 
to couple to the v = 17 bound state of the 0g (P3/2) poten- 
tial which is at an energy A = —6.85 x 10 4 r from dissocia- 
tion. For a lattice with wavelength 785 nm and transverse and 
parallel confinements V x = 38. 5£^, V" = 22Er, the result is 
F v ~5x 10~ 7 . Given this confinement, the on-site interaction 
using a s = 5.6 nm is U = 3.574kHz. Choosing an experi- 
mentally reasonable atomic Rabi frequency Q.a — 25T, where 
r = 2% x 6.065MHz, we find K w 0.1317. Here the atomic 
scattering due to the catalysis laser is saT w 6.7 x 10~ 8 r per 
atom and the off resonant light shift is |V C .| = 15.5U. 

By way of example we define a one dimensional regis- 
ter ^,with 501 atoms that resides inside a lattice filled with 
N = 551 atoms. An external magnetic trapping frequency 
of (07- = 2k x 8Hz ensures that the last occupied well has 
an energy e((N — l)/2) = 0.9U. We note that the proba- 
bility for tunneling of holes from the edges is negligible as 
J I («8) = 0.28. In practice it is not important to know the ex- 
act number of atoms in the lattice as long as the trap strength 
is chosen such that, given the uncertainty in the number of 
atoms, the constraint e((N — l)/2) < U is always satisfied. 
These parameters fix the ratio U/J — 500 and the measure- 
ment strength is therefore K./2y/nJ w 1.5. 

When the environment is monitored, for instance by 
looking for photon scattering from the bound molecular 
state, the evolution of ground states can be modeled us- 
ing quantum trajectories. For our simulation, the ground 
state wavefunction |v|/(f)) = cr(t)\T) +Y.j,± c s ± ( t )\Sf) i s U P _ 
dated according to the non-Hermitian Hamiltonian H =Hj — 
ijM/2Hj,± ^ quantum trajectory corresponding 

to a null measurement result converges to the target state \T) 
and freezes the state there as demonstrated in Fig. |2](for a sim- 
ilar effect with ions see Q)- The preparation time scale is 
tprep = 1/k. The probability that the preparation fails, i.e. the 
system decays from the excited states is pf a u = 1 — pr^O). 

Real experiments will be constrained to finite detector ef- 
ficiencies T|. For T) = 0, corresponding to nonselective mea- 
surement, the system dynamics evolve according to Eq. [9] If 
we represent the dynamics of the system as a pseudo two state 
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FIG. 2: Target state fidelity, F(t) = pr,r(0. forU/J = 500 and n = 
500 atoms. The main plot shows a quantum trajectory for a state 
that begins in the ground state of the Bose-Hubbard Hamiltonian and 
converges to the unit filled state \T) by a null result on a selective 
measurement. The time scale for saturating the target state for the 
parameters here is t sa , ~ K _1 = 7.7 /U. The inset shows the evolution 
for three detector efficiencies (r| = 1.0,0.9,0.8) over a longer time 
scale. Also shown is the oscillatory dynamics with frequency U if 
the measurement is turned off after the target state is reached. 



system \T) and \S), with an average energy splitting U + \V C \ 
the equations of motion for the pseudo Bloch vector are: 



u 

v 
w 
x 



(U + \V c \)v-xu _ 
-kv - (U+ \V c \)u - y/2nJw 
-x(x + w) + A\f2nJv 
-k(x + w), 



(11) 



where u = Re[ps.7-], v = Im[ps j-], w = ps.s — Pt.t, and the 
decreasing norm is x = Tr[p]. After a period 1 /|V C |> the coher- 
ences approach steady state, and the target state population, 
assuming pr,r(0) ~ 1 and k/2 v / «7 > 1, is 



PfAt) = PT,r(0)e 



-(&nJ 2 Kl)/((U+\V c \] 



(12) 



This behavior exhibits the continuous quantum Zeno effect 
in the limit K./2y/nJ — > oo. Even for moderately strong mea- 
surements, if the initial state is close to the target state then 



the deviation time is long compared to \/U. For finite effi- 
ciencies, we can express the approximate fidelity to be in the 
target state. Assuming a null measurement result, for times 

t > t prep ,l/\V c \, it is F(tv) = Pr,r(f) = n + -i\)Pt,tW- 
In practice, high detection efficiencies may be obtainable by 
applying a second photo ionizing laser on resonance with the 
molecular bound state and monitoring the emission of ions. 

If the catalysis field is turned off after the target state is 
reached, the system will freely evolve according to Hbh- The 
general time dependent state in the restricted basis is: 



\V{t))=c T {t)\T) + Y j e 
;",± 



-iE(Sf)t 



f(t)\Sf 



(13) 



Assuming the initial state |W(0)) = \T), the amplitudes 
ct,c s ± evolve according to the Schrodinger equation, 

i^W*)) = HBH^it)). Using first order perturbation the- 
ory the time dependent fidelity to be in the target state 
F(t) = K^l 1 ^?))! 2 can be estimated under the assumption that 
S/U < 1. For times t < 1/7 we find 

F(t) = l-8(7/f/) 2 («-cos((7f)(l+sin(8(n-l)f)/sin(50). 

(14) 

This solution compares well with exact numerical simula- 
tions for a small number of atoms (N,M < 8) and simulations 
for larger N = M using a restricted basis set of dimension 
iV(iV — 1) + 1 consisting of states with at most one site with 
two atoms. Note that the time averaged deviation from perfect 
fidelity 1 — (F(t)) = 8«(7 /U) 2 is twice as bad as the deviation 
if the system were left in the ground state Because the 

Bose-Hubbard Hamiltonian is intrinsic to trapped bosons in 
an optical lattice, other dispersive techniques for initiali zing 
a register in a lattice such as Raman side-band cooling 1 15], 
and phase space compression 1 16], if not corrected, will suffer 
from the same loss of fidelity as described by Eq. El 

In summary, we have shown that efforts to prepare a reg- 
ister of atomic qubits in an optical lattice suffer from errors 
inherent in the underlying many body dynamics. We have 
introduced a protocol that addresses this issue to make the 
MI transition a robust mechanism for initialization. While we 
have presented the idea in the context of one dimensional dy- 
namics, the method is also applicable to higher dimensions, 
which is the relevant regime for scalability. We appreciate 
helpful discussions with Almut Beige, Paul Julienne, and Eite 
Tiesinga. This research was supported in part by ARDA/NS A. 
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